Fokker-Planck equation for transport of wave packets in nonlinear disordered media 
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Starting from first principles, we formulate a theory of wave packet propagation in a nonlinear, 
disordered medium of any dimension, through the derivation of a Fokker-Planck transport equation. 
Our theory is based on a diagrammatic expansion of the wave packet's density, and is supported by 
a heuristic picture that involves a Boltzmann equation with an effective, external potential. Our 
approach also confirms results obtained in previous work for two-dimensional, nonlinear disordered 
media. 
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I. INTRODUCTION 



When it propagates in a disordered environment, a particle (or more generally a wave) usuallyexperiences a classical 
random walk, with ballistic motion interspersed with a large number of scattering processes The overall motion 
is then described by a diffusion equation for the ensemble-averaged particle's density. In some cases however, this 
simple picture breaks down due to constructive interference that build up between multiple scattering paths, leading 
to alteration of the diffusive description. Weak localization, and in the extreme case Anderson localization, which 
corresponds to a complete halt of transport, are examples of these interference phenomena (T], 0. Another effect 
that may strongly affect transport is the possible presence of nonlinearities. In the optical context for instance, the 
nonlinearity can be obtained by increasing the intensity of the light propagating in the disordered medium, making use 
of the sensitivity of the refractive index to the electric field (the Kerr effect) . Recent works exploited this phenomenon 
to study the interplay between disorder and nonlinearity in systems of optical fibers 0, Q • 

Another promising system for studying quantum transport in disordered environments and in the presence of 
nonlinearity are ultracold atomic gases. In these systems, the disordered medium is produced by means of an optical 
random potential, and the nonlinearity originates from interactions between atoms. In this context, the particular 
problem of a Bose-Einstein condensate prepared in an optical trap and then released in a random potential has 
■ sparked considerable interest in the last few years [5[. The major advantage of this type of experiment is that 
the time and spatial evolution of the macroscopic wave function of the condensate can be observed by imaging the 
i atomic density profile. Furthermore, experiments involving condensates offer unprecedented control of the physical 
parameters driving the disorder and the nonlinearity. 
00 , Until now, much experimental effort has been devoted to the study of Bose-Einstein condensates in one-dimensional 
random potentials, with strong interest in the phenomenon of Anderson localization, the role of atomic interactions, 
C*~) , and the interplay between these effects [f|. In recent works, Anderson localization of a condensate was experimentally 
demonstrated, in a situation where interactions were presumably negligible @)H|]- Although one might expect that 
interactions play no role after a long period of expansion of the condensate, during which it becomes increasingly 
dilute, recent theoretical works suggested that, even at long times, Anderson localization could be destroyed by 
interactions, leading to a subdiffusive regime [9l4ll|. 

In higher dimensions, our knowledge of the behavior of Bose-Einstein condensates in the presence of disorder is still 
in its fledgling stage. Among the handful of experiments that have been carried out, one can cite recent beautiful 
works involving ultracold atoms, in which the Anderson transition was demonstrated and analyzed in detail, through 
the experimental realization of a quantum-chaotic system, the quasiperiodic kicked rotor, known to be equivalent to 
the three-dimensional disordered system [l2|, [l3| . Regarding the quantum transport of Bose-Einstein condensates in 
random potentials, the common feature between the two- and three-dimensional cases is that unlike in one dimension 
(ID), a broad diffusion regime exists, provided the disorder is not too strong. Here, "diffusion regime" means that most 
atoms of the condensate propagate by diffusion. In three dimensions (3D), it was shown that this situation is achieved 
if the quantity fir/H is large [14]], where jj, is the chemical potential of the condensate and r the scattering time, i.e., 
the average time between two successive scattering processes. In 3D, both the diffusion [l5[ and the localization [l6[ 
regimes have been studied theoretically, as well as the cross-over between them [bH ]. From the experimental point 
of view, the diffusion of a cloud of thermal atoms in a random potential was recently studied [ijj ■ Yet, in two and 
three dimensions, whether it concerns a regime where localization effects dominate or a situation where transport is 
driven by diffusion, the role of atomic interactions on condensate expansion remains largely unknown. Furthermore, 
to our knowledge no analytical description of this type of problem is available. In a recent theoretical work however, 
a kinetic equation describing the propagation of a wave packet in the infinite two-dimensional disordered medium was 
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proposed, based on the introduction of an effective nonlinear potential [l8| . From this equation, the authors of Ref. 
[18| were able to infer the qualitative behavior of the wave packet in the presence of nonlinearity, and, in particular, 
put forward the phenomenon of "locked" explosion, corresponding to a global diffusive behavior, but with a rapid 
explosion at early stages of the expansion, located in the central part of the wave packet. 

This paper aims at formulating a consistent theory able to describe wave packet propagation in a nonlinear, weakly 
disordered medium. For this purpose, we derive a Fokker-Planck equation describing transport of such a wave packet, 
starting from first principles. This equation holds in any dimension, and, in 2D, reduces to the kinetic equation 
introduced in Ref. jl8| . Our approach is based on a diagrammatic treatment of the wave packet's density in the 
presence of nonlinearity. Such a formalism was initially introduced to describe the stationary transport of waves 
in nonlinear media either formed by nonlinear point scatterers at random positions (l9l [20| or by linear scatterers 
embedded in a homogeneous nonlinear medium (2lj . In particular, it was successfully applied to the description of 
the coherent backscattering effect in the presence of nonlinearities [22J. In this paper, we extend this approach to 
time dependent diffusion processes, and, in addition, validate the diagrammatic approach by a heuristic treatment 
that consists in solving the Boltzmann equation in the presence of an effective, nonlinear potential, in the spirit of 
[l8j . We work in the framework of quantum particles, but our results can be extended to classical wave scattering 
described by the nonlinear Helmholtz equation with a disordered potential [21| . 

The paper is organized as follows. In Sec. II, we introduce the physical quantities relevant to describe transport of a 
wave packet in a nonlinear disordered medium. We introduce the Fokker-Planck equation and establish the connection 
with the work in Ref. [18| . The two main approximations used throughout the paper, the diffusion approximation 
and the limit of weak disorder, arc discussed, and the validity of the latter in a nonlinear medium is analyzed in Sec. 
III. In Sec. IV, we recall the core elements of the derivation of the standard diffusion equation describing wave packet 
dynamics in a linear disordered medium. This derivation is extended to the nonlinear case in Sees. lyland fVTl leading 
to the Fokker-Planck equation. Finally, in Sec. IVII1 we show how this equation can be straightforwardly obtained 
from a heuristic approach resorting to the Boltzmann equation in the presence of an external, effective potential 
describing the nonlinearity. Technical details of calculations are collected in the three appendices. 



II. THEORETICAL FRAMEWORK 



Before examining the nonlinear case, we first recall some basic facts concerning the propagation of a wave packet 
in the infinite, linear disordered medium of dimension d (d = 1, 2 or 3). This wave packet can represent, for instance, 
a collection of particles of mass m moving in a disordered environment. At a given point r inside the medium and at 
time t, the wave function of this wave packet obeys the time-dependent Schrodinger equation: 
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2m 



V(r) 



(i) 



In Eq. ([T]), and in the rest of the paper, we set H = 1. The disordered medium is described by the random potential 
V(r), that we assume to obey the white-noise Gaussian statistics 



V(r)V(r') = 7<J(r - r'), 



(2) 



where the overbar denotes averaging over realizations of disorder and 7 = 1/(2ttut), with r the scattering time and 
v = (d/2)(m/27r) <i / 2 e d ' /2 ~ 1 /r(l + d/2) the density of states per unit volume. By Fourier transforming Eq. ([T]) with 
respect to time we obtain 



2m 



V(r) 



&(r) = 0, 



(3) 



where i/'e( r ) = / dt ip(r, t)e z<it . In this paper, we will be interested in the average density at energy e, defined as 

n e (r,t) = J ^-n t (r,u)e-^\ (4) 



where n e (r,cj) = tl) e+u / 2 (r)ijj*_ u , 2 (r). Physically, n e (r, t) can be interpreted as the probability density of finding a 
particle with energy e in the vicinity of point r, at time t. Even in the absence of nonlinear effects, the evaluation 
of this object is, in general, a difficult task that requires to average over the product of the sum of all partial waves 
reaching the observation point r. The analysis is, however, greatly simplified if one restricts to the so-called diffusion 
approximation and to the limit of weak disorder [23j . In the diffusion approximation, one considers only large spatial 
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scales |r| 3> I (assuming the wave packet is initially located at point r = 0), where i is the mean free path, and slow 
dynamics u <C e, 1/r. The limit of weak disorder corresponds to the condition er 3> 1 and describes complete phase 
randomization of the partial fields, such that no interference survives the disorder average [23]. "Phase coherent" 
effects such as weak or strong (Anderson) localization are thus neglected in this picture. Under these assumptions, 
the average density obeys the diffusion equation [23| 

(d t - D t V 2 )n c (v, t) = 5(t)n e (v, t = 0). (5) 

The right-hand side of this equation is the "source term" , described as a short pulse whose magnitude is given by 
the wave packet's density at the initial time. The main lines of the derivation of Eq. (|5|) will be recalled in Sec. IIVI 
D e = v 2 t Jd is the diffusion coefficient at energy e in dimension d, expressed as a function of the scattering time and 
the velocity v — \/2e/m. From the solution n e (r,t) of Eq. ([5]), the total density of the wave packet is obtained after 
integration over energies: 

f d 

n(r,i) = |</>(r,t)| 2 = J ^n e (r,t). (6) 

Let us now introduce nonlinearity in the medium. If this nonlinearity is "not too strong" , the wave function of the 
wave packet obeys the Gross-Pitaevskii equation 



-^V 2 + V{v)+g\^{v 1 t)\ 2 



^(r,f). (7) 



The term g\ip\ 2 plays the role of a nonlinear potential, repulsive if the parameter g is positive, and attractive if g 
is negative. In the context of ultracold atomic gases, the Gross-Pitaevskii equation controls the time-evolution of a 
Bose-Einstein condensate, described in terms of a macroscopic wave function ip, in the presence of a random potential 
V(r). The nonlinearity then originates from interactions between atoms [25j | . 

The presence of the nonlinear term in Eq. (JT)) substantially complicates the analysis of the problem. In this paper, 
we show that in a nonlinear, disordered medium, the diffusion equation ([5]) for the average density n e (r,t) must be 
replaced by a Fokker-Planck equation that reads 

[dt - VD e _ e V + d t 9d e + V^gF] n e (r, t) = <5(*)n e _ e (r, t = 0), (8) 

where we have introduced the "nonlinear potential" 8 — 8(r,t) — 2gn(r,t) and its associated "force" F = — V0. Eq. 
([5]) is valid in any dimension, and thereby generalizes the kinetic equation analyzed in [18j . which only applies to the 
two-dimensional case (for the sake of clarity we have kept the same notations as in [l8|). 

Although Eq. ([8]) formally describes the propagation of a wave packet in a disordered, nonlinear environment, 
its analysis reveals that the problem is analogous to the one of a collection of independent particles, moving in a 
disordered region and subject to the external force F = —"V9. Each particle has total energy e and kinetic energy 
e — 8. The parameter r] t -e = (d/2— l)/x e _g depends on the dimensionality d and on the quantity /i e _e = Z? £ _g/(e — 8), 
which can be seen as the mobility of a particle with kinetic energy e — 8. The main consequence of the presence of 
the potential 8 is that the diffusion coefficient, which depends on the kinetic energy e — 8, now acquires a position 
and time dependence. The third term in the left-hand side of Eq. © originates from the time dependence of the 
potential 8, and the last term corresponds to a drift effect. For d = 2, we have r\ t — and one recovers the kinetic 
equation analyzed in [l8j]. Although Eq. ([S} is rather complicated, we can gain more physical insight by performing 
the substitution eo = e — 8, and integrating the resulting equation over eo- We then obtain the following simple 
continuity equation for the total density: 

d t n(r, t) + V ■ J(r, t) = 5{t)n{r, t = 0), (9) 

where the total current J is defined as 

/d 
l-D to Vn eo (r,t) +M £o Fn eo (r,t)] . (10) 

The continuity equation guarantees the conservation of the total density n(r, t). The formula (| 10[) has a clear physical 
interpretation: the total current is given by the sum of a diffusion current (present in the linear case), and a drift 
current, induced by the force F. This drift current is proportional to the mobility ji eo — D eo /eQ. Again, Eq. (jlOp 
holds in any dimension. In the particular two-dimensional case, D e cx e, such that Eqs. © and (fTU)) can be simplified 

as 



d t n(r,t) = -V 2 
m 



J ^ e ° rle o( r ' t ) +9n 2 {r,t) 



(11) 
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for t > 0. From the particular simple form of Eq. (|TT|) . it follows that the mean radius of the wave packet grows as 
t 1 ! 2 in 2D, similarly as in the linear case. This phenomenon, named "locked explosion" in Ref. [l8|, is a consequence 
of the fact that the mobility fx eo is independent of eo in 2D. In ID or 3D, we expect different physical scenarios because 
the mobility acquires an energy dependence. This is however still an open question that we leave for later work. 

The Fokker-Planck equation ([8} is a promising tool to study transport in nonlinear disordered media, and therefore 
calls for microscopic justification. Derivation of this equation from first principles and the formulation of a consistent 
theory of dynamical transport in a nonlinear, disordered medium, are the two purposes of the present contribution. 
In the following, we will focus on the three-dimensional infinite medium, and explain how results are modified in ID 
and 2D. 



III. LIMIT OF WEAK DISORDER IN THE PRESENCE OF NONLINEARITY 

In the microscopic approach developed hereafter, we will make use of the limit of weak disorder er 1 discussed 
in Sec. [TTl assuming that this condition is not modified in the presence of nonlinearity. This assumption amounts 
to neglect the effect of the nonlinearity at the scale of the scattering time, and leads to a restriction for the values 
taken by the average density. To explicit this restriction, we express that the scattering time tnl associated with the 
nonlinear potential gn(r, t) in Eq. ([7]) should be much larger than the scattering time r associated with the disordered 
potential |26l |27|. To calculate tnl, we use the short-range "Ci" part of the spatial correlation (5n e (r, t)(5n e (r', t) of 
density fluctuations [23]. The condition tnl S> t then leads to: 

T 5 n(r,;)« v^- (12) 

Condition ([T2l will be assumed to hold from here on. Note that to obtain it, we neglected the long-range "C2"part 
of the correlation function of density fluctuations (23[. This approximation is justified in our case, since one readily 
shows that at long times, C2 oc (l/i^ 1 )/ (kt) d ~ 1 in dimension d [l4j]. The time decay of C2 then allows us to neglect 
its contribution to tnl in 2D and 3D. In ID, this assumption is questionable because Ci <~ C\. This is related to 
the fact that in ID, localization effects take place at very short time scales t ~ r, which, strictly speaking, limits the 
validity of the diffusive approach. We will come back to this point in Sees. IVII and [Villi 

IV. DIAGRAMMATIC APPROACH IN THE LINEAR CASE 

In this section we recall the derivation of the usual diffusion equation ([5]) in a linear, disordered medium. The 
microscopic description of transport is based on a Feynmann representation of the density. More precisely, in the 
limit of weak disorder, this description is a perturbation theory which consists in finding the leading diagrams in 
l/(er) contributing to n e (r, oS) = ip ei (r)-0* 2 (r), where ei = e + lo/2 and e 2 = e — w/2. These diagrams are formed by 
combining a scattering sequence associated with ip ei and a scattering sequence associated with ip* . The key point 
is that in the limit of weak disorder, the only diagram that survives the disorder average is the so-called "ladder" 
diagram, for which ip ei and ip* 2 follow exactly the same scattering sequence, thereby introducing no phase difference 
between ip ei and t/j* . The ladder diagram is depicted in Fig. [TJi. It connects the scattering sequence associated with 
tp ei (solid line) with the one associated with ip* 2 (dashed line). Inner dotted lines symbolize the correlation function 
of disorder @ between a scattering event experienced by ip ei and one experienced by tp* 2 . The ladder diagram is 
generated by the "building block" depicted in Fig. Qi>, whose mathematical representation is the Bethe-Salpeter 
equation [231 ]: 

V> ei (r)^(r)=? 6l (r)£ 2 (r) + 7 J dV G €l (r', r)^ (r', r)^ (r')< (r'), (13) 

where G e (r',r) is the average "amplitude Green's function", which is the Green's function of the Schrodinger equation 
([T]), averaged over disorder. In Fig. \Tjp, G ei and G £2 are depicted by the solid and dashed lines, respectively. The 
Bethe-Salpeter is an integral equation for the density. In the diffusion approximation and in the limit of weak disorder, 
it can be simplified to yield the diffusion equation. For this purpose, it is convenient to work in Fourier space. We 
thus introduce n € (q,w) = J d 3 r ip ei (r)ip* 2 (r)e~ lq ' r and S e (q,ui) = J d 3 r ip ei (r)i/; e2 (r)e~ l ' i ' r . A calculation detailed in 
Appendix [Al then gives 



n e (q,u) ~ 5 e (q,w) + (1 + iuir - D e rq 2 ) rc e (q, w) 



(14) 
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FIG. 1: (a) Ladder diagram, connecting the scattering path associated with by ip ei (solid line) to the scattering path associated 
with tjj* (dashed line), r is the observation point. Arrows indicate the direction of propagation. Inner dotted lines symbolize 
the correlation function of disorder given by Eq. (b) Building block generating the ladder diagram, and leading to 

the Bethe-Salpeter equation (|13[) . The solid and dashed lines depict Green's functions G ei and G E2 , respectively. The last 
scattering process of the ladder diagram occurs at point r'. (c), (d), (e) and (f) Building blocks generating the contributions 

^> n (r)^iJ a (r)|cj ip ei { r )i , S 2 ( T )\<ii ^<si ( r )V'?2 ( r )l<= an d ^nWV'ejWIf to tne density, respectively. The black square represents the 
nonlinear parameter g. Other symbols have the same meaning as in diagrams (a) and (b). 



where D e — v 2 r/d is the diffusion coefficient at energy e. Taking the inverse Fourier transform of Eq. (fl4| with 
respect to cj and q, we obtain 

(d t -D e V 2 )n e (r,t) = ^^-, (15) 

T 

which is nothing but a diffusion equation with source term 5 , e (r,t)/r. A calculation detailed in Appendix iBl leads to 

S e (v,t)=TS(t)n e (r,t = 0), (16) 

where, in the limit of weak disorder and within the diffusion approximation, the wave packet's density n e (r,t = 0) at 
the initial time is given by (see Appendix [Bj 

n e (r,t = 0)~ jj 27 ^ QMQ + q/ 2 )<TO - q/2)e-" q - r - (17) 

In Eq. (|17l) . A(e, Q) = — ImG c (Q)/7r is the "spectral function", which is the density probability of finding a particle 
with energy e among all particles with momentum Q (l6j . In the limit of weak disorder er ^> 1, the spectral function 
becomes very close to the free-space expression A(e, Q) ~ 5(e — eq), where eg = Q 2 / (2m) [l^, |T3] . </>(k) is the Fourier 
transform of the wave function cf>(r) of the wave packet at the initial time. The integral over q accounts for the 
possible finite extent of the wave packet. For the particular case of a "point source", <f>(r) oc S(r) and consequently 
S e (r,t) oc d(r)S(t). In this case, n e reduces to the Green's function of the diffusion equation, sometimes called 
"probability of quantum diffusion" in the literature [23j . 



V. DIAGRAMMATIC APPROACH IN THE NONLINEAR CASE 



A. Nonlinear diagrams 



In the previous section, we saw that in the linear case, only the ladder diagram contributes to the density in the 
limit of weak disorder. In the nonlinear case, each of the two scattering paths building the ladder diagram may be 
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"perturbed" by the nonlinear potential g\ip\ 2 in Eq. (J7J, at any point of the scattering sequence. The density at some 
point r is then obtained by summing an infinite series of diagrams of the type of the one depicted in Fig. [5] In this 
figure, every ladder diagram originates from a point of the wave packet at the initial time, and is connected to one of 
the two scattering paths forming another ladder diagram, via two average amplitude Green's functions. From here on, 
we will refer to such a connection as a "nonlinear scattering process" . In Fig. [2 the black square at each connection 
symbolizes the interaction parameter g. More details about this construction can be found in }20l |. 
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FIG. 2: A typical series of diagrams contributing to the density in the nonlinear case. All ladder diagrams constituting this 
series originate from a point of the initial wave packet, and are connected to another ladder diagram via two average amplitude 
Green's functions. This connection is referred to as a nonlinear scattering process in the main text. Only the central diagram 
is connected to the observation point r. 



Since an infinite number of diagrammatic sequences like the one of Fig. [2] must be accounted for in order to calculate 
the density, writing down an integral equation similar to the Bethe-Salpeter equation of the linear case seems to be a 
difficult task. This is however possible if we notice that all possible diagrammatic series are composed of the only five 
"building blocks" depicted in Figs. [Tb, Ht, HH, QJ andQJ. The density is obtained by summing the integral equations 
associated with each of these building blocks. This strategy was initially developed in the particular case of a finite 
nonlinear disordered medium under plane wave excitation, i.e., corresponding to cj = (or equivalently £i — £2) t 20]. 
In the present case, an additional difficulty proceeds from the energy exchange arising at each nonlinear scattering 
process. This exchange has to be examined carefully when writing down the integral equations for diagrams [IJd, [TJ;, 
QJl, [TJb and QJ. It can be readily understood from the Gross-Pitaevskii equation which, in Fourier space, can be 



rewritten as 



14] 



2m 



The nonlinear term in the left-hand side of Eq. (1181) has a straightforward interpretation: when propagating in the 
disordered medium, the partial field V>£-£i+e 2 1S affected by the density ^ £l ?/i* 2 , from which an amount of energy ei — £2 
is transferred. In general, this energy exchange may affect either the field tp or its complex conjugate tp*. This is 
manifested by the existence of the two building blocks in Figs. [T}1 and QJ, respectively obtained from the building 
blocks in Figs. QJ: andQJ, after the substitutions G O G , t\ -O- £2 and £3 O £4. 

Before proceeding further, two important comments are in order. First, keeping all the building blocks in Figs. 
QJ), QJ, QJi, QJ: and QJ is crucial in order to account properly for nonlinear effects, while preserving conservation of 
n(r,t). Second, although the present approach sums infinitely many diagrams, one has to keep in mind that it remains 
perturbative since we do not consider situations where more than one nonlinear scattering process occurs in a row. 
This means that we implicitly assume that the parameter rgn(r, t) is small, such that condition (TT21 is automatically 
fulfilled in the limit of weak disorder er> 1. As we will discuss in Sec. IVIII however, the Fokker-Planck equation (j8]) 
resulting from the diagrammatic approach is valid beyond this assumption. 



B. Calculation of building blocks 



We now examine the integral equations corresponding to the building blocks in Figs. [TJ: and [TJi. The one in Fig. 
[Th generates a series of diagrams for which the scattering sequence associated with ip ei emerges at point r from a 
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nonlinear scattering process, and the scattering sequence of from a linear scattering process. From here on we 
denote by tp ei (r)ip* (r)| c the contribution to the density associated with this series of diagrams. It obeys (see Fig. 
Hfc): 

de% de 



^ 1 (r)^ 2 (r)| c = 2 7ff // dVd 3 n // ^^G^(i^,r)G 61 _ w+M (i^,r 1 )(? ei (r 1 ,r)^ es (r 1 )<(r 1 )x^ 1 _^(^*Jr'). 



Note that the densities that appear inside the integrals, and evaluated at points r' and ri, must correspond to the total 
density in order to effectively sum all diagrams. The prefactor 2 accounts for two possibilities to group the incoming 
Green's functions converging at point ri into pairs (see Fig. [TJ:). Similarly, the building block of Fig. [TJi generates a 
series of diagrams for which the scattering sequence associated with ip ei emerges at point r from a linear scattering 
process, and the scattering sequence associated with tp* from a nonlinear scattering process. The associated series of 

diagrams V^i( r )*fe( r )|d obeys the same equation as ^> ei (r)^?| 2 (r)| c , but with the substitutions G G , t\ «-> e 2 and 
£3 f> £4. In the diffusion approximation and in the limit of weak disorder, the integral equations for the diagrams [TJ; 
and[l]i can be simplified. This calculation is lengthy, and is reported to Appendix [C] for clarity. It can be conveniently 
performed in Fourier space, where it leads to 



. ffdEdfl f d^ 

AgT JJ 2^2^ J T^Wqi^XCq-q^-O) 

ff dEdn f d 3 qi n . . 
+ gT J J 2^2^ J 7^ fi ^( q i^)9 £ n £ (q- qi , W -0), 



2 



2k£ 2k£ 



^(q -q qi> 



(20) 



where k = y / 2me, n e (q,w) = J d 3 r i/j ei (r)ip* 2 (r)e _lq ' r , and we recall that e = (ei + £2)/2 and ui = t\ — e 2 . 

The building blocks in Figs. QJ and [If are calculated in the same way. For instance, the diagram QJ yields the 
integral equation 

V> ei (rX(r)| e = 2 7 2 5 j j j d 3 v' a9v x afv 2 ff (r', r)G £l _ £3+£4 (r', r 2 )G £l _ £3+£4 (r 2 , n)G £l (r 1; r 2 )G £l (r 2 , r) x 



X^e 3 (rX(ri) x ^ £l _ £3+£4 (r)V £ * 2 (r'), (21) 

where we used Eq. ^ to account for the additional impurity line. Once again, the integral equation for the building 
block in Fig. QJ follows from Eq. (|2"Tj) after the substitutions G -O- G , ei -O- ei and €3 ■<-» £4. Summation of the 
building blocks e and f leads to (see Appendix [C]) 



n e (r,u;)|e + n e (r,u;)|f = 4 5 t J J 7^-7^ J 4^n s (q 1 ,0)n e (q- q 1: w- fi) x 



1 + 2z(o/-Q)t r , , q?/ 

— 2^ — + 2^ ( - q +q ' qi -y) 



(22) 



In order to complete the diagrammatic treatment, one finally needs to include the building block in Fig. [TJd, which 
describes a usual linear scattering process. This type of process may of course also occur in a nonlinear medium (see 
Fig. [2]). When this single building block is taken into account, we saw in Sec. IIVI that it generates a ladder diagram 
described by the Bethe-Salpeter equation Eq. (fT"3|) . In the nonlinear case, it generates a series of diagrams for which 
both scattering sequences associated with tp ei and tp* emerge at point r from a linear scattering process. This series 
of diagrams contributes of an amount ip ei (r)'ip* (r)\ h to the total density, and obeys 



^e 1 (r)V* 2 (r)|b = V ei (r)^ 2 (r)+7 J d 3 r' G £l (r', r)G £l (r', r)^ £l (r')V £ * 2 (r'). (23) 

The simplification of Eq. ([23]) follows exactly the same lines as for the Bethe-Salpeter equation (see Appendix ■ 
By analogy with Eq. ([14")) , we have 

n 6 (q,w)| b = S' £ (q,w) + (l + ia;r-D £ Tq 2 )n £ (q,a;), (24) 

where we have defined the source term 5 £ (q, to) — J d 3 r ip ei (r)ip e2 (r)e _lq r . The prime symbol in this definition signals 
that in the nonlinear case, the source term has not the same value as in the linear case. The expression of 5* £ (q, w) 
will be given in the next section. 
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C. Combination of nonlinear diagrams 

we have to close Eqs. I®, ([22]), and ([24 



In order to complete the derivation of the Fokkcr-Planck equation 
This is achieved by writing 

n e (q,uj) = n e (q,w)| b + ra £ (q,w)| c + n e (q,u)\ d + n e (q,u))\ e + n e (q,U))\{. 

Using Eq. d20j) , ([22]) and ([24]), we rewrite Eq. ([25]) as 

-S' e (q,Lj) = (iujT - D £ rq 2 )n e (q, lj) + 

+4.gr 



(25) 



dE dn r d 3 qi 

27T 27T ./ (27T)' 3 



i(w-Q)r £ 2 q 2 
2fc^ 2fc^Y 



/•/• f d 3 q, n . _ . 



(26) 



We now divide Eq. (|26[) by r, and take the inverse Fourier transform with respect to q and d. We obtain 



-L> e VX(r,i) + 

2gf 
3k£ 



2gr f dE , ,' 



„•> f dE , ' 



n e (r,t)+2g 



d t n £ (r,t) - 

^ f dE , . 
/ 2^ UE ^ r,t ' 



d e n e (r,t). 



(27) 



From Eq. ([2"T1) . we see the emergence of the "potential" 9(r,t) = 2g J(dE/2Tr)riE('r,t) = 2gn(r,t). Since 2rgn/(k£) ~ 
2rgn/(eT) <C 1, we can divide both sides of Eq. (|27p by 1 + 2rgn/(k£) and expand the result for small nonlinearity. 
Using £ 2 /(3k£) = Dj(2e) in 3D, we can rewrite Eq. ([27]) as 



= _ ^ e v 2 + d t 6(r,t)d e ] n € (r,t) + ^ [0(r, t)W(r, i) - u e (r, i)V 2 #(r, t)] 



(28) 



In order to complete the proof, we have also to evaluate the source term which appears in the left-hand side of Eq. 
([2"5]l . Although this term describes a propagation free of linear scattering processes, it is likely to be affected by the 
nonlinearity [20|. The result (TTT)) obtained in the Sec. IIVI should therefore be reconsidered here. The calculation of 
S"(r, t) is achieved by summing a series of diagrams of the type in Figs, lc and Id, but where the correlation functions 
ipijj* of the field at points r' and ri are replaced by the correlation function (jxff of the initial wave function of the 
wave packet, and where an arbitrary number of nonlinear scattering processes may occur between the points r' and 
r. This calculation is lengthy and is reported to Appendix [B] for clarity. It yields: 



5'(e,r,t) = r<y(t)ne_j(r,t = 0), 



(29) 



which is the same expression as in the linear case but with the total energy e replaced by the kinetic energy e — 0. 
The density n e _g(r,t = 0) of the wave packet at the initial time appearing in Eq. ([2T)| is still given by Eq. (TTT)) . but 
the spectral function A must now be evaluated at the energy e — 9. We now combine Eqs. ([28]) and ([29l) . and slightly 
rearrange the terms in the right-hand side. This yields 



$(t)n e _ 9 (r, t = 0) = [d t + d t 6d e ] n e (r, t) - V 



D, 



-27) Vn e (r,t) 



§(V0)n e (r,i) 



(30) 



If we notice that D e - D e 6/(2e) ~ £> e „ e and Dj(2e) = (d/2 - l)D e /e ee r/ e in 3D, we see that Eq. (0Qj is nothing 
but the Fokker-Planck equation ([5]) in the limit gTn(r,t) <C 1, which completes the microscopic derivation. 



VI. DIAGRAMMATIC TREATMENT IN ID AND 2D 



In this section we briefly discuss how the above derivation is modified in ID and 2D. The building blocks to consider 
remain evidently the same as in the three-dimensional case (see Figs. [IK [It, Gib BJ and [if), but their calculation leads 
to different expressions. The two-dimensional case was already sketched in [181 ]: the last two terms in the right-hand 
side in Eq. ([30]) are then replaced by 



D, 



[V0(r, t) Vn e (r, t) + 0(r, t) V V(r, t)] 



(31) 
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which yields the Fokker-Planck equation © because in 2D, r\ e = 0, and -D £ _g = D e — D e 9/e. Finally, for completeness 
we also analyze the one-dimensional case, although it is less relevant due to the existence of strong localization effects 
that take place at time scales of the order of the scattering time [Hj]- In ID, the last two terms in the right-hand side 
in Eq. ([50)1 are replaced by: 



e 



2V9(r,t)Vn e (r,t) + ~0(r,t)VV(r,t) + ^n e (r,t)V 2 9(r,t) 



(32) 



In ID, ?7 e = —D e /(2e), and D e _g ~ D t — 3D e 9/(2e) in the limit rgn(r,t) <§; 1, such that Eq. (|32|) leads once again to 
the Fokker-Planck equation. 



VII. HEURISTIC DERIVATION OF THE FOKKER-PLANCK EQUATION 

In the microscopic derivation presented in the previous sections, the "nonlinear potential" #(r, t) — 2gn(r, t) emerges 
naturally. Note that 9 differs by a factor 2 from the term gn(r, t) added to the potential V in the Gross-Pitaevskii 
equation [[?]), Apart from that, it turns out that the complicated problem of wave packet propagation in a nonlinear 
disordered medium can be replaced by the more simple one of a collection of independent, classical particles subject 
to an external force, and moving in a disordered environment. The purpose of this section is to confirm this picture, 
and to take advantage of it to present a very simple way of deriving the Fokker-Plank equation l[5]). 

For a system of independent particles subject to the "external" force F = — V6>, the transport properties are 
contained in the one-particle density function /(r, p, t), which gives the probability density of finding a particle at 
position r with momentum p and kinetic energy eo = p 2 /(2m). Since we are considering a disordered medium with 
point-like scatterers, the time evolution of / is governed by the Boltzmann equation 



dt + - ■ V + F ■ V p 

m 



/(riPlt) = M«. (33) 



Eq. (f3"3"|) expresses the fact that the rate of change in / in phase-space (r, p) is due to the scattering from impurities. 
The density n £o (r,t) introduced in the previous sections is related to / through the simple relation n eo (r,t) = 
is f(r,p, t), with v the density of states per unit volume. We stress at this point that eo does designates the kinetic 
energy, which should not be confused with the total energy e of Sees. [V] and EH f(r,p,t) = J dil p /{Aii)f{r, p, t) is 
the average of /(r, p, t) over all directions Q p of p. In order to find a transport equation for f(r,p, t), we make use of 
the diffusion approximation, which here consists in expanding the density function into spherical harmonics, keeping 



only the lowest two harmonics [28[: 



/(r,p,t) ~ f(r,p,t)+dn p ■ h(r,p,t), (34) 

where f\(r,p, t) = J df2p/(47r)ri p /(r, p, t), and d — 1,2,3. Substituting this expansion into Eq. (|3"3")l . we obtain an 
equation, which, after integration over angles, yields 



d t f(r,p,t) + 



rad rad 60 ( d) p 



fi(r,p,t) = 0, (35) 



and, after multiplication by tt p and integration over angles, 



fi (r,p, t) = ( V + Fd C0 ) /(r, p, t) . (36) 



Finally, combining Eqs. (1351) and (1361) . using the relation n eo (r, t) = v f(r,p, t) and writing the initial value n eo (r, t = 0) 
explicitly in form of a source term in the right-hand side, we obtain 

d t n eo (r,t) - [V + Fey D eo [V + Fd eo - ( V£ jD eo )F] n eo (r,t) = 5(t)n t0 (r,t = 0), (37) 

where i] eo = (d/2 — l)D eo /eo. Eq. (|37|) was obtained in [18| in the particular two-dimensional case, for which the 
parameter rj t0 — 0. Eq. (|37|) is the generalization of this result to any dimension. The Fokker-Planck equation 
(jHJ follows straightforwardly after the transformation eo e ~ 9 (with e the total energy), in agreement with the 
diagrammatic treatment. Note that in this section we only made use of two approximations. First, we resorted to the 
diffusion approximation, which allowed us to use the expansion (1341) . Second, when writing the Boltzmann equation 
(1331) , we assumed that the scattering time was unaffected by the nonlinearity, which is valid as long as condition (|12[) 
is fulfilled. Unlike in the previous sections however, we did not use the condition Tgn(r,t) <C 1. This means that Eq. 
(13"T1) . as well as the Fokker-Planck ((SJ), are not restricted to this limit. 



10 



VIII. CONCLUDING REMARKS 

In this paper, we have formulated a transport theory of wave packet propagation in a nonlinear, weakly disordered 
medium, through the microscopic derivation of the Fokker-Planck equation ((5J) . The latter describes transport in media 
of any dimension, and thereby generalizes the equation introduced in [181 ] . Our derivation is based on a diagrammatic 
formalism, under the assumption of weak disorder, and in the framework of the diffusion approximation. Within this 
approach, we have identified all the diagrams required to ensure conservation of the total density n(r, t), and we have 
derived the Fokker-Planck equation in the limit of small nonlinearity (rgn <C 1). As compared to the linear case, the 
main difficulty of the nonlinear diagrammatic procedure presented in this paper is that nonlinear terms always show 
up in powers of the parameter rgn/(er). This can be seen, for instance, in Eq. (1281) . where the nonlinear terms dt0d e , 
(-D e /e)#V 2 and (D c /e)V 2 # in the right-hand side are all of order Tgn/(er). The perturbation theory has therefore to 
be carried out carefully, keeping terms up to the first order in rg , n/(er). Note that the fact that nonlinear corrections 
to the diffusive transport depend on the single parameter rgn/ '(er) means that it is not easily possible to decouple the 
nonlinearity from the strength of the disorder within a perturbative approach. This effect was already pointed out in 
the case of the three-dimensional expansion of a Bose-Einstein condensate in a disordered medium, in the presence of 



very weak interactions [14 



We have also shown that the physical problem of wave packet expansion in a random potential can be conveniently 
described by a heuristic approach, in which the wave packet is replaced by a collection of independent, classical particles 
subject to an external, effective potential. This approach, based on the Boltzmann transport equation, allowed us 
to recover the Fokker-Planck equation, thereby validating the diagrammatic procedure. Unlike the latter however, 
the derivation from the Boltzmann equation is non-perturbative, and only relies on the diffusion approximation and 
the fact that the scattering time is not affected by the nonlinearity. This means that the Fokker-Planck equation (|SJ) 
remains valid for possibly large values of rgn, on condition that rgn does not exceed \fer. 

Let us conclude our study by a discussion on the practical validity of the Fokker-Planck equation in the time domain. 
In deriving Eq. (J5J , we have assumed that no coherent effects of localization type take place in the disordered medium. 
This assumption limits the range of validity of the description at long times, when the mean radius of the wave packet 
reaches the localization length £ £0 , where e is the typical kinetic energy of the wave packet. For the particular case 
of a Bose-Einstein condensate expanding in a random potential after release from an optical trap, eo is of the order 
of the chemical potential \i [lE| . In ID and 2D, we can estimate the time t\ oc after which the approach breaks down 
due to localization effects, by requiring that (r 2 ) ~ £ 2 , where (r 2 ) ~ Dfj_t\ oc . In ID, this yields ti oc ~ r, which means 
that the diffusive approach breaks down at short times, and thus limits the relevance of the present work in this case. 
In 2D, t\ oc ~ r exp(/rr) 3> r, such that a broad diffusion regime exists where the Fokker-Planck equation applies. In 
3D, the situation is slightly more complicated since there exists a critical energy, separating diffusive and localized 
states. It was shown recently that in this case, the diffusive approach, and consequently the Fokker-Planck equation 



is valid, at distance r from the origin, as long as time does not exceed ti oc ~ r(r/£) 8 / 3 (er) 2 ^> r 14|. 
The transport equation © is a flexible theoretical tool, and we think it could be used as a basis for future work on 
time-dependent problems in disordered and nonlinear media, for instance in the context of Bose-Einstein condensates 
expanding in random potentials. Although strictly speaking we have only addressed transport in the infinite medium, 
Eq. ([SJ holds in finite media as well, provided it is supplemented by adequate boundary conditions [23 |. Since the 
nonlinear effects only weakly affect transport at the scale of one mean free path, we expect the boundary conditions 
used in the linear problem to remain approximately valid in nonlinear media. 
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Appendix A 

In this appendix we present a derivation of Eq. (TT4"1) . starting from the Bethe-Salpeter equation (|13p . For this 
purpose, it is convenient to work in Fourier space. The Fourier transform of the second term in the right-hand side 
of Eq. (TT3)) is given by 

d 3 k — _ * 

gG ei (q)G ea (k + qK(q,w), (Al) 



(2tt 
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where n e (q, u>) = J d 3 r ip ei (r)ip* 2 (r)e lqr , with e — + £2)/2 and w = ei — £2- We have also introduced the Fourier 
transform of the average amplitude Green's function |23[ 

e- e q + i/(2t) 

with e q — q 2 /(2m). In the diffusion regime, |q^| <C 1, and we can expand the Green's function at energy £2 according 
to 

C; (k + q) ~ (k) + [(v • q) + e q ]G: 2 (k) 2 + (v • q)^ (k) 3 , (A3) 

where v = q/m. Inserting this expansion into Eq. (jAll) . and performing integrations over k, we obtain 

f d 3 k — — * v 2 t 2 
7 J ^3G £l (q)G £2 (k + qK(q, W ) ~ (1 + twr - -^-q 2 H(q,w). (A4) 

This result, combined with Eq. (TTS]), yields Eq. (TH| . with D e = u 2 r/3. In ID and 2D, the derivation follows exactly 
the same lines, the only difference being the replacement of the diffusion coefficient by D e = v 2 r/d. Note that in 
deriving Eq. (|A4I) . we neglected all terms of order l/(er) or higher (limit of weak disorder), and kept only terms 
linear in ujt and q 2 £ 2 (diffusion approximation). 

Appendix B 

In this appendix we calculate the "source terms" S e (r,t) and S' e (r,t) in the diffusion and Fokker-Planck equations 
(IT51) and (j2"8)) . respectively. Both are defined as 

SP(r,t) = J ^ e+u/2 (C u/2 (r)e- fat . (Bl) 

To begin with, we calculate S e , i.e. the source term in the absence of nonlinearity. To do so, it is convenient to relate 
the wave function tp at point r to the wave function <j> of the wave packet at the initial time through 

lk(r) = |d 3 r s G e (r s ,r)0(r s ). (B2) 

Substituting Eq. (|B2|) into Eq. ([HTjl. we obtain 

Se(T,t) = J ^~ tUlt JJ rf 3 r Sl d 3 r S2 G £W2 (r Sl ,r)G;_ w/2 (r S2 ,r)0(r Sl )^(r S2 ). (B3) 

The diagrammatic representation of Eq. (|B3[) is shown in Fig. Introducing the Fourier transforms G(q) = 

J (( 3 rG(r)e" lq ' r of the average amplitude Green's functions appearing in Eq. (|B3[) . we obtain 

S e (r,t) = J ^ lUt JJ (2^3 ff 3 Ge +u /2 (Q + q/2) G* e _ u/2 (Q - q/2) <j> (Q + q/2) (/>* (Q - q/2) . (B4) 

We now perform the integral over w. Using |q| <C |Q|, uj <C e (diffusion approximation), and er > 1 (limit of weak 
disorder, we find: 

J ^^Ge+u.MQ + q/2)G £ *_ w/2 (Q - q/2) ~ 2-KT8(t)A(e, Q), (B5) 
where we have introduced the spectral function A(e, Q) = — ImG c (Q)/7r. This yields: 

5(6, r, t) = r*(t) I f ^ J^^e, Q)^ (Q + q/ 2 ) ^* (Q - q/ 2 ) e ~ iqr - (B6) 

Eq. (|B6|) is the source term that appears in the diffusion equation (fT"5|) . It can be readily related to the wave packet's 
density at the initial time if we note that n e (r, t = 0) = J(dw/27r)?/; £+w / 2 (r)7/>*_ w / 2 (r) — J(c?o;/27r)^ e+tJ / 2 (r)'0 e _ w / 2 (r). 
Following exactly the same reasoning as above, we find 

n e (r,t = 0) = || ^^2 7 rA( e ,Q)0(Q + q/2)^(Q-q/2) e -^, (B7) 
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FIG. 3: (a) Diagram describing the source term in the absence of nonlinearity. r 31 and r S2 are two points of the wave packet 
at initial time, (b) and (c) Leading-order diagrams contributing to the source term in the presence of nonlinearity, in the limit 
rgn -C 1. (d) Typical diagram contributing to the source term, allowing to account for large values of the nonlinearity. 



which yields immediately Eq. (1161) of the main text. 

We now consider the source term S' e , which is affected by the nonlinearity. In comparison with the linear case 
where only the diagram in Fig. [5^, is relevant, now also the nonlinear diagrams in Figs. andJSJ; contribute. For the 
sake of consistency with Sec. |Vj we assume that Tgn(r, t) <C 1, namely we neglect diagrams for which more than one 
nonlinear scattering event occurs in a row. The final result that we obtain is however valid beyond this limit, as we 
discuss later on. In the nonlinear case, the source term is given by 

S' £ (r, t) = S e (r, t) + S® ( r , t) + S^(r, t), (B8) 

where the two terms with the superscripts (b) and (c) refer to the diagrams in Figs. [3]d and[3j;, respectively. Within 
the diffusion approximation where lot, (63 — e±)T <C 1, is for instance given by 

S^(r,t) ~ 2gJ ^e- iUt J[ f d^dh^r, JJ ^^^.^(n.^lr^r) x (B9) 
x <j>(r sl )4>*(r S2 ) x ^(ri)^* (r 2 ). 



We now approximate V , e 3 ( r i)V , * 4 ( r i) — ^esMVCM) an d perform the integrals over e 3 , 64, w and 1*1, by using the 
following result, valid in the limit of weak disorder: 



f d d nG e (r Sl , n)G 6 (n, r) ~ -it'- ^G e (r Sl , r). (BIO) 

Eq. (IB9|) then becomes 

5( fe )(r,t) = -S(t) [[ d 3 r si d 3 r S2 2lT ^ ~ rsi 1 n(r ,t = 0)G t (r ai , r)^ (r, 2 , r)0* (r Sl )0(r S2 ). (Bll) 



The calculation of S^ c ' is identical. It leads to the same result as Eq. (IB11[) . but with |r — r si | replaced by — |r — r s 
Then, with Eq. JB3|, Eq. (JB8]) yields 



S' e (r,t) = S(t) // d 3 r Sl rf 3 r S2 exp(-Ai|r-r Sl |)exp(Ai|r-r S2 |)G e (r Sl ,r)G;(r S2 ,r)0(r Sl )0*(r S2 ), (B12) 



where A = 2gTn{v, t = 0)/£. Note that to obtain Eq. (|B12|) . we made use of the approximation 1 — 2igTn\r — r Sl \ /£ + 
2igTn\v — r S2 \/t ~ exp(— 2igTn\v — r Sl \ j I + 2igTn\r — r S2 1). This approximation, valid in the limit rgn <C 1, turns 
out to be the exact result obtained when one assumes that an arbitrary number of nonlinear scattering events may 
affect the scattering paths associated with G and G* , i.e., when one evaluates all diagrams of the type of the one in 
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Fig. [3ji. Eq. (|B12[) therefore holds also for large values of rgn. By introducing the Fourier transform of the Green's 
functions and the exponentials that appear in Eq. (|B12|) . we rewrite S' e as: 

SKr.t) - i(i )////|i||^0B e (Q- qi+q /2)B:(Q + qi - q /2)x (B13) 

xF_(q 1 )F+(q 2 )#(Q + q/2)<(,'(Q - q/2)e-"", 

where -F=p(q) = f d 3 r exp(=FiAr — iq • r) = (27r) 3 <5(<7 ± A)(5(cos0)/<7 2 , with 6* the angle between q and the z-axis. 
Performing integrals over q x and q 2 , and using the diffusion approximation |q| <C |Q|, we obtain 

S' £ (r,t)=6(t) jj ^_^|G e (g + A)| 2 0(Q + q/2)r(Q-q/2) e ^ r , (B14) 

where Q = |Q|. By definition [see Eq. (|A2I) ]. we have 

|G £ (Q + A)| 2 = -2rImG e _ £Q+A (Q) =27rr^(e-e Q+A ,Q), (B15) 

where we have used the definition of the spectral function A in the last equality. The last step consists in writing 
A(e — eq + x,Q) — A(e — £q — 2Ae/fc, Q) = A[e — €q — 2gn(r,t = 0),Q], which is a good approximation as long as 
Tgn(r, t) <C er. This finally leads to the same expression (|B6I) as in the linear case but with the energy e replaced by 
e — 8 in the spectral function, and completes the proof of Eq. ([29)) . 



Appendix C 

In this appendix we derive Eqs. (|20[) and (|22|) . starting from Eqs. (fT9|) and (|2Tj) . respectively. To do so, we introduce 



the notation n(q, ei, £2) = n e (q, u) = J d 3 rip ei (r)i/>* (r)e * qr , with e = (ex + £2)/2 and w = ei — e 2 . 
Consider first n e (r, w)| c . Taking the Fourier transform of Eq. ((T5)) with respect to r, we obtain 



n e (r,cj)| c = 27.9 



d 3 q a d 3 k f f de 3 de^ 
(2tt) 3 (2tt) 3 



2tt 2tt 



G* 2 ( k - q)G , ei - e3+£4 (k-q 1 )G ei (k)n(q 1 ,e3,e 4 )n(q- q^ei - e 3 + e 4 ,e 2 ) 

( C1 ) 

In the diffusion approximation, |q£|, \q_i£\, |(ei — £2)7"! "C 1 and | (£3 — &i)t| <C 1. The integral over k can then be 
performed by expanding the Green's functions G e2 and G £l _ C3 + e4 in powers of q and q 1 , by means of Eq. (| A3|) . This 
leads to 

d 3 k 



, 2 N3 G e 2 ( k -q) G e 1 -e3+6 4 (k-qi)G ei (k) = 
-» + 2(a,-n)r- 1 + ^- n ) T + W-q. qi ; 





i^qfl 




3 



(C2) 



where f2 = e 3 — £4. As we pointed out in Sec. |Vj n e (r,w)|d can be straightforwardly deduced from n e (r, cj)|d via the 
substitutions G -f-> G , ei O £2 and £3 f> £4. With this procedure, we can write 



n e (r,w)| c + ra e (r,w)| d = 
d 3 qi 



2g T 



-2gr 



(2tt) 3 
j - 
d 3 qi 



27^ 2vr' 



n £; (q 1 ,0)n e _n/2(q- qi,w - fi) x 



1 + - fi)r 



2fc^ 



(2tt) 3 77 2vr 2tt 
i- 2(w- fi)r - 



1 + i(u> — fi)r .0, 2 . 

■ + £ 2 (q 2 -q- qi ) 



/J__ \_^qf 



2fc^ 



V 2M 



(C3) 



where we changed the integrations over e 3 and 64 into integrations over E = (e 3 + 64) /2 and = e 3 — £4, such 
that n(q 1 ,e 3 ,e 4 ) = n E (qx,fi), n(q - q^ex - e 3 + £4,^2) = n e _n/2(q- qx7 w - and n(q - q 1 ,e 1 ,e 2 + £3 - £4) = 



14 



n e+n/2{Q — — The next step consists in expanding the densities ?if±o/2 around e, making use of fl <C e 
(diffusion approximation). We thus have 



(C4) 



Substituting this expansion into Eq. (IC3|) . we obtain 

n e (r,w)| c + ra e (r, ui)\ d = 

. ff dEdtl f rf 3 qi , n . . . 

dEdn f d 3 q x 



i + i(ui-n)T i 2 , , 

2fc^ + 2W (q - q ' qi) 



2 ^ // TT-TT- 



2?r 2tt 7 (2tt 



r fin £ ;(q 1 , f2)<9 e n e (q- q^w - 0), 



(C5) 



which is Eq. (|2"0)) of the main text. Note that in deriving Eq. (IC5|) . we kept only first-order terms in wt, Or, q 2 ^ 2 
and q 2 ^ 2 (diffusion approximation), and in rgn/(k£) ~ Tgn/{er). Terms of higher-order have been neglected. 
We now examine n e (r,w)| e . The Fourier transform of Eq. (I21[) with respect to r reads 



n e (r,w)| c = 27 2 5 



,/'q, ,/'k ,/'k' /■^^^^_^_^^_^ )Bai _ <i+ai(k ._ qi)5ai(k)B<i(k / )x 



(2tt) 3 (2tt) 3 (2tt) 3 
X«(qi, £3, e 4 )n(q - q 1( ei - e 3 + e 4 , e 2 ) 



(C6) 



As compared to Eq. (jCll) . an additional integral over k' comes into play. At the leading order in l/(k£) ~ l/(er), 
this integral reads 



7 



d 3 k' 

(2tt) 3 



2H' 



(C7) 



The remaining integral over k is the same as the one in Eq. (|C1[) . Keeping only first-order terms in Tgn/(k£), lot, 
fir, q 2 £ 2 and q 2 ^ 2 , we have 



7 

IT 

2k£ 



d 3 k d 3 k' — * — — — — 

(2tt) 3 {2-Kf Gt2 ^ ~ ^ G ^-^+^ k ~ <h) G ei-<*+u ( k ' - qi)G £l (k)G £l (k') 

-z + 2(w - n)r + J?Mq 2 - q • q x + 



(C8) 



Using the same substitution as above to evaluate n e (r,w)|f, we obtain 

n e (r,w)| e + ?i e (r,w)|f = 
rf 3 qi 



2gr 
2gr 



(2tt) 3 

rf 3 qr 
(2tt) 3 



dE dVL . . i r . , , q, 

y o-^Efai! ")"e-n/2(q - qi, w - O)— |-« + 2(cj - Q)t + ri(q z - q • q : + y ) 

dEdfl , . , — i r . . 9 q, ' 

y j-^lqii ^K+n/2(q - Qi, w - fi) — |z - 2(w - Q)r - r«(q' i - q • q x + y ) 



(C9) 



We now expand the densities within the integrals, using Eq. (|C4[) . However, this time we can restrict ourselves to the 
zero-order term in the expansion (|CJ4|) , since keeping the term proportional to Q would lead to terms of second-order 
in lot, fir and q 2 £ 2 , which are negligible in the diffusion approximation. We thus write n e ±o/2 — n e , and obtain 



n 6 (r,w)| e + n e (r,u)\f 



*9T // 7: t; 



di? dO /• d 3 q 1 



2?r 2tt J (2?r) 3 
l + 2i(w-Q)r £ 2 



n jE (q 1 ,0)n e (q- q x ,w - 0) x 



2fcg 



^(-q^q-q^f) 



(CIO) 



which is Eq. 



of the main text. 
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